The roughness of stylolites: 
Implications of 3D high resolution topography measurements 
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Abstract 

Stylolites are natural pressure-dissolution surfaces in sedimentary rocks. We present 3D high 
resolution measurements at laboratory scales of their complex roughness. The topography is shown 
to be described by a self-affine scaling invariance. At large scales, the Hurst exponent is £i ~ 
0.5 and very different from that at small scales where £2 ~ 1.2. A cross-over length scale at 
around L c = 1 mm is well characterized. Measurements are consistent with a Langevin equation 
that describes the growth of a stylolitic interface as a competition between stabilizing long range 
elastic interactions at large scales or local surface tension effects at small scales and a destabilizing 
quenched material disorder. 

PACS numbers: 83.80.Ab, 62.20.Mk, 81.40.Np 
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Stylolites are geological patterns that are very common in polished limestones, a material 
largely used to construct floors and walls of buildings and monuments. They are observed 
as thin irregular interfaces that look like printed lines on rock cuts, which is responsible 
for their name. They are roughly planar structures that are typically perpendicular to the 
geological load (i.e. lithostatic pressure or tectonic maximum compressive stress). These 
rock-rock interfaces are formed at shallow depths in the Earth's crust during deformation of 
sedimentary rocks and result from a combination of stress-induced dissolution and precipita- 
tion processes! They a r e found in m any sedimentary rocks such as .unestones, sandstones 



or evaporites [2| and exist on a very large range of scales, from micro-meters to meters. 

Despite their abundance, stylolites are, as mentioned by Gal et al. [jj, "among the least 
well- explained of all pressure-solution phenomena". First they are complex 3D structures 
that are often only described from 2D cross-sections since they are generally partially sealed 
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. Second, they develop in various geological contexts which lead to very different ge- 
ometries. Third they are sometimes transformed because of processes like diagenesis and 
metamorphism that develop after their initiation. 

In this Letter we show the first 3D high resolution topography measurements of natural 
stylolite interfaces that could be fully opened. We characterized the scaling invariance, 
namely self-affinity, of the morphology and show the presence of a cross-over length scale. 
We also propose a model of the stylolites roughening. It is based on a Langevin equation 
that accounts for stress-induced dissolution in a quenched disorder. 

The roughness measurements have been performed on three independent stylolite inter- 
faces included in very fine-grained limestone samples from Burgundy area, Vercors, and Jura 
mountains in France (Fig. [IJ. The samples have been collected in newly open quarries, thus 
preserved from late breakage and chemical erosion. The opening procedure was possible for 
these samples because of the accumulation of undissolved minerals like clays that formed a 
weak layer along the stylolite interface. The concentration of these minerals provides an es- 
timate of the cumulative strain by dissolution the sample underwent . As shown in Fig. d 
peaks along the interface are randomly distributed in space and of various sizes (up to one 
centimeter). Large peak magnitudes and local high slopes along the topography makes the 
roughness measurement difficult and challenging. 

We used two different profilo meters to sample the stylolite roughness. First, with a 
mechanical profilometer [f], ~\ we extracted four profiles of 1030 points each with a horizontal 
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FIG. 1: Picture of a stylolite surface (S12A) in a limestone from Vercors Mountains. Magnitude 
of the peaks are typically of the order of 6 mm. 
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FIG. 2: A ID profile obtained by a mechanical profilometer (1030 data points - Ax = 30/um) 
along a stylolite surface. 



step of Ax = 30 /xm. The mechanical profilometer measures the surface height from the 
contact of a needle onto the surface. The radius of the needle tip is 25 /im. The vertical 
resolution is 3 /im over the available range of 5 cm. One of these profiles is shown in Fig. [2j 
We compare the mechanical measurements to an optical profiling 8]. This technique is based 
on a laser triangulation of the surface without any contact with the surface. The laser beam 
is 30 /im wide. Horizontal steps between measurement points were Ax = Ay = 7 to 50 fim 
with a vertical resolution of 2 /im. The main advantage of this technique comes from the 
high acquisition speed that can be performed compared to the mechanical profilometer, since 
there is no vertical move and on-flight measurements are possible. However, a successful 
comparison with mechanical measurements is necessary to ensure that optical fluctuations 



are height fluctuations and not material property fluctuations. Three independent samples 
have been measured at very high resolution: one side of a stylolite from Jura mountains 
(Sjura) with a resolution 600 x 600, one side of a stylolite from Burgundy area (S15) with a 
resolution 8200x4100 and two opposite surfaces of the same stylolite from Vercors mountains 
shown in Fig. □ with a resolution 2400 x 1400 for S12A and 8200 x 4100 for S12B. 

We analyzed the height distribution in terms of self-affinity [9| which states that the 
surface remains statistically unchanged for the transform: Ax — > A Ax, Ay — > A Ay, 
Az — > A*° Az, where A can take any real value. The exp onent ( is the so-called Hurst 
exponent. A ID Average Wavelet Coefficient technique |li| has been used. For a self-affine 
profile, the wavelet spectrum behaves as a power law with a slope 1/2 + (, and provides an 
estimate of the Hurst exponent (. The spectra clearly exhibit two regimes (Fig. (3J) ■ At large 
length scales, a power law behavior is observed with a slope of 1 in the log- log plot, which 
is consistent with a Hurst exponent of (j = 0.5. At small length scales, a second power law 
behavior is observed with a larger slope (1.7) in agreement with a Hurst exponent £2 = 1.2. 
The crossover length scale is sharp and defines a characteristic length scale which is slightly 
different for the three surfaces, L c »s 1 mm. L c is several orders of magnitude larger than 
the grain size and significantly larger than experimental cutoffs. This spectral behavior is 
observed for both mechanical and optical measurements. 

We checked that another analysis technique, namely the Fourier power spectrum, was 
providing very consistent results. Fig. 0] shows averaged ID spectra of profiles extracted 
from the surface Sjura. A self-affine property of the profiles leads to a power-law behavior 
of the power spectrum as P(k) oc fc _1_2 ^J3]. Moreover, average spectra of profiles taken 
along perpendicular directions provide very consistent results (Fig. Ej). Isotropy of scaling 
invariance is confirmed by the circular symmetry of the 2D power spectrum of the surface. 

The second part of the letter is devoted to a modeling of the stylolite roughening. The 
aim is to understand the origin of the self-affine behaviors and of the characteristic length 
L c . We propose to consider a simple model as a paradigm for an interface growth in an 
heterogeneous medium like natural rocks. This model aims at providing a framework for 
future modeling of stylolite formation. 

We consider the following geometry: The stylolite interface is assumed to be initiated 
along the boundary between geological beds. Accordingly it can be approximated as the 
boundary of a quasi-flat and very elongated fluid pore. The trapped fluid is assumed to 
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FIG. 3: Averaged wavelet spectra of topographic profiles extracted from four optical maps of 
stylolite surfaces. Spectra have been normalized to superimpose for large L on the spectrum of 
S12A. 



form a film and to be at lithostatic pressure. The solid, where this pore is embedded, is 
supposed to undergo an average stress cr° = <j zz zz + a® x (xx + yy), where z refers to the 
direction normal to average stylolite direction and x and y refer to directions along the 
average stylolite direction. Since stylolites are on average normal to the largest principal 
stress direction, a s = \<J° ZZ \ — |cr° x | > 0. 

Possible solid contacts with the mirror surface on the other side of the fluid film are 
neglected, considering that such contact points concentrate strain when they occur, and 
induce faster dissolution of these contacts, thus leading to an essentially lubricated contact 
zone between neighboring grains: for simplicity, we neglect interactions between the mirror 
surfaces and assume that the front morphology is to first order dominated by a dissolution 
process between a fluid film and a single elastic solid. 

Assuming a free surface profile z(x,t), the normal n to the interface pointing toward 
the solid is, in the limit of small relief, h = z — (d x z)x, where we assume plane strain 
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FIG. 4: Fourier power spectra of ID topographic profiles oriented along two perpendicular direc- 
tions (X and Y) and of the full 2D surface. The latter was radially integrated to be compared to 
the ID power spectra. Inset shows a gray map of the 2D power spectrum. A mirroring technique 
has been used to reduce non periodic edge effects. 



perturbations. The stress state in the solid is expressed as a = a + a , where mechanical 
equilibrium between solid and fluid requires that cr 1 • (— n) = — a s (d x z)x. This stress state 
results from a surface distribution of tangential force -a s (d x z)x applied on the quasi-planar 
boundary of the solid by the fluid, so that using Green's elastostatic function and 
integrating along the y— direction, at the surface, a xz = a\ x = a s (d x z) and a xx = a yy = 
a s {2v /tt) f dy(d y z(y)) / (x — y), all other components being null. 

For small reliefs ( | |/||co| I ^ 1) an d to leading order the elastic free energy u e = 
[(1 + v)OijOij — vakk&u}/ can be approximated as u e = u° e + u\ where from the above, 

< = apl/E (1) 
u\ = - (3( Po a s /E) J dy(d y z)/(x - y) (2) 

with an average solid pressure po = — (2a xx +a^ z )/3, and two dimensionless positive constants 
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a = [9(1 - 2v) + 2(1 + v)a 2 s /pl}/12 and (5 = u(l - 2v)/n, where E is an effective Young's 
modulus, and v the Poisson coefficient. u° e is the elastic energy from the global tectonic 
loading and u\ is its local perturbation that results from the interface topography. 

The chemical potential difference at the solid/fluid interface that can potentially desta- 
bilize the interface can be written as jl2|: 

Afj, = Q(u e + 7«) (3) 

where u e is the elastic energy per unit volume in the solid, 7 is the surface energy, k the 
curvature, and Q a molar volume. We have assumed that gravity effects are negligible. 
We have also assumed that the matrix of the solid, i.e. an assembly of initial sedimentary 
particles, is sufficiently porous during stylolites initiation to have a bulk diffusion within 
the material. This assumption is supported by rock thin section observations under an 
optical microscope [l3| . If a bulk diffusion holds in the fluid surrounding the stylolite, the 
evolution of the interface is directly related to the chemical potential: v n = mA/j, where 
v n is the normal dissolution velocity and m is the mobility Q|. We also neglected the 
chemical potential evolution within the film since we only aim at describing the initiation 
of the process under drained conditions. 

This homogeneous description thus predicts, for small reliefs, d t z = Vq + mVL{u\ + , yd xx z), 
with vq = mVtu®. Surface tension is a stabilizing term, but it is important to note that 
the elastic interaction term, ul, is also stabilizing in the present context. For the present 
situation, stylolites are perpendicular to the maximum principal stress, and will subsequently 
be assumed horizontal: a s > 0. Considering an elementary departure from a flat interface, 
such as a fluid intrusion in the solid, i.e. a bump with a maximum in x, such as d y z > 
for y < x, and d y z < for y > x, ul is negative in x and reduces the dissolution speed in 
the bump at x. Accordingly, since the problem is linear, elastic interactions are stabilizing 
for any corrugations of the interface. For vertical stylolites, the picture is opposite (a s < 0) 
and elastic interactions are destabilizing leading to a lateral expansion of the stylolite. 

The homogeneous picture predicts the propagation of a planar dissolution interface driven 
by the average elastic energy u° e) with an average speed estimated as v ~ 8 • 1CT 6 m/year 
where we used m = kfl/(RT), with a dissolution rate k w 10 -4 mol/m 2 /s, fi ~ 4 • 10~ 5 
m 3 /mol for calcite, R is the universal gas constant, T w 300 K, a ~ 0.5, E sa 8 • 10 10 Pa for 
limestones, a characteristic stress estimated as po ~ 25 MPa, corresponding to a rock at 1 
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km depth. 

To understand the dynamic roughening of stylolites, it is essential to capture the effect 
of heterogeneities of relevant material properties in the solid, namely u, E, m and 7. We 
assume the relative variation (5E/E and others) of these properties to be small, and to 
correspond to independent random variables associated to each constitutive grain of the 
rock, which are typically £=10 /im sized. At early stages of the process where d x z < 1, we 
define the dimensionless surface position with respect to the average plane z' = (z — Vot)/£ 
and the dimensionless space and time variables x' = xjl and t' = t/r where r = £ 2 / { r yVtm) 
to obtain, to leading order in relative fluctuations and typical slopes, for the roughening 
interface speed: 

d t ,z'(x', = d x , x ,z> -j-f dy'^- l + rj(x', z'{x')) (4) 

L J x — y 

where L* = ^E/(l3p a s ) and 77 = [a£p /((3L*a s )} ■ [(SE/E) + (Sm/m) — {5a/ a)}. In this 
Langevin equation with quenched noise, the destabilizing random term is balanced by the 
restoring surface tension term at scales below L*, and by the restoring elastic interactions 
at scales above L*. We propose that this critical scale L* corresponds to the measured 
crossover length L c . For typical limestones, 7 = 0.27 J/m 2 for a water-calcite surface and 
v « 0.25, so that (3 ~ 0.04 and L* « 0.9 mm, consistently with the above measured. The 
other characteristic quantities of interest are r ~ 0.2 year and the characteristic amplitude 
of the dimensionless noise rj is p ~ aipo/ (j3\*a s ) ~ 0.2. 

For the Laplacian regime (L L*) and the mechanical regime (L ^> L*), only one of the 
two restoring terms in Eq. Q dominates, and these two independent regimes have already 
been studied. Indeed, the Laplacian regime is nothing else than the Edwards Wilkinson 

n 

(EW) problem [14| in a quenched noise. In this case the interface is self-affme with an 
exponent £ 2 ~ 1.2 [l^. In the mechanical regime, Eq. (jlj) is analogous to the quasi-static 
propagation of an elastic line or a mode I fracture front in a disordered material, and the 
Hurst exponent is Ci ~ 0.4 for a kernel similar to Eq. (j2J) 16, 17, "Lsl \v\. 

The roughening amplitude can be obtained by considering the EW equation with 
quenched noise regime: the characteristic width of the surface measured at scale L scales 
as w(L)/£ w p(L/£) ( ' 2 at saturation, obtained from a flat interface after a saturation time 
r s (L) such that t s /t « (L/£)&/ s , with a dynamic exponent 5 w 0.8 [15]. With L » 1 
mm, I fa lOyum and (2 ~ 1-2, this scaling law predicts up to a constant of order unity the 
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saturation width at cross over scale w(L*) rs 0.5 mm and the time to saturation as t s m 200 
years. This length scale corresponds to the measured one (Fig. |2J, and the short saturation 
time implies that observed stylolites have achieved their saturation width over geological 
time scales. That the width amplitude is also correctly predicted in the mechanical regime 
could be checked directly, but is granted by the fact that it is correctly predicted in the 
Laplacian regime, as well as the crossover scale, which determines entirely the prefactor of 
the scaling law w(L) in the L > L* regime. In principle, determining L* and w(L*) could 
give two independent constraints on both p and a s , which could allow to determine both 
the pressure and differential stress prevailing during the formation of a particular stylolite. 
However, given the amount of approximations in the involved constants, the only way to 
test this effect on the cross-over wavelength would be to measure stylolites formed in various 
geological conditions and study the effect of depth and orientation to the main stress. 

In conclusion, we presented a quantitative description of stylolite interfaces. The ex- 
perimental measurements are 3D high resolution descriptions of the topography of natural 
stylolites. We show that the surfaces are self-afnne but with two regimes. At small scales, 
the Hurst exponent is unexpectedly high, (2 = 1-2, and consistent with a Laplacian regime. 
At large scales, the stylolites morphology is controlled by long range elastic stress redistri- 
butions. In this case the roughening is important with a low Hurst exponent (1 = 0.5. The 
two regimes are separated by a crossover characteristic length L c , also predicted by a model 
based on the description of a stress-induced dissolution, where restoring surface tension ef- 
fects and elastic interactions compete with a quenched noise. It is important for geological 
implications to note that L c is very sensitive to the average stress pq. Indeed, a measurement 
of L c from roughness profiling could provide an estimate of the stress magnitude during the 
stylolite growth, that is, in the past. Accordingly stylolites could be considered as stress 
fossils. 
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Perfettini, P. Meakin, and E. Merino for fruitful discussions and two anonymous reviewers 
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